On the effective velocity created by a point vortex in two-dimensional hydrodynamics 
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We complete previous investigations on the statistics of ve- 
locity fluctuations arising from a random distribution of point 
vortices in two-dimensional hydrodynamics. We show that, 
on a statistical sense, the velocity created by a point vortex 
is shielded by cooperative effects on a distance A ~ n~^^^, 
the inter-vortex separation. For i? A, the "effective" veloc- 
ity decays as instead of the ordinary law recovered 
for R <^ A. These results are similar to those obtained by 
Agekyan [Sov. Astron. 5 (1962) 809] in his investigations 
on the fluctuations of the gravitational field. They give fur- 
ther support to our previous observation that the statistics of 
velocity fluctuations are (marginally) dominated by the con- 
tribution of the nearest neighbor. 

PACS number(s): 47.10.-fg, 47.27.-i, 47.32.-y, 02.50.-r 

Recently, several papers [|l|- Q have derived the form 
of the velocity probability distribution function (p.d.f) 
created by a random distribution of point vortices in two- 
dimensional hydrodynamics. Technically, the velocity V 
occuring at the center O of a circular domain of radius 
R is the sum of the velocities $i {i = 1, TV) produced 
by the N vortices 
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where denotes the position of the ith point vortex rel- 
ative to the point under consideration and, by definition, 
r±i is the vector rotated by -|-7r/2 (we have assumed for 
simplicity that all the vortices have the same circulation 
7). Therefore, if the vortices are randomly distributed 
over the entire domain, the problem at hand amounts 
to finding the distribution of a sum of random variables. 
The difficulty is that the variance of the individual ve- 
locities ($^) diverges logarithmically so that the Central 
Limit Theorem is not directly applicable. In early works, 
Jimenez , Min et al. and Weiss et al. |^ determined 
the velocity p.d.f by applying a generalization of the Cen- 
tral Limit Theorem due to Feller |^ and Ibragimov & 
Linnik |ll|] . More recently, this problem was reconsidered 
by Chavanis & Sire and Chavanis |^ who extended 
the methods introduced by Chandrasekhar & von Neu- 
mann in their investigations on the fluctuations of 
the gravitational field. In the stellar context, the limit 
distribution is a particular Levy law, called the Holtz- 
mark distribution. In the vortex context, the velocity 
p.d.f. is given by 
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{V<Vcr^u{N)), (2) 
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where 



{V>Vc„t(N)), 
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VcUN)^ \^-^\nNj \n'/\\nN). (4) 

This distribution is obtained in the "thermodynamical" 
limit N,R ^ -\-oo with n — finite (we have assumed, 
for simplicity, that the vortices are uniformly distributed 
on average; see Refs. [§,||j9j for a generalization of these 
results to the case of an inhomogeneous medium). In 
fact, this limit is not well-defined and the velocity p.d.f. 
is "polluted" by logarithmic corrections. Therefore, in 
Eqs. (||)(||)(|4|) we must consider that -f 00, but not 

IniV. This is appropriate to physical situations in which 
the typical number of point vortices does not exceed 10*. 
Since this distribution is intermediate between Gaussian 
and Levy laws we proposed to call it a "marginal Gaus- 
sian distribution" [Q. In the (formal) limit In ^ +00, 
the tail is rejected to infinity and the limit distribution is 
Gaussian [|10|jll[ . However, the convergence is very slow 
so that, in physical situations, the algebraic tail is clearly 
visible §1). 

From the distribution (||) (^) (^ we can easily calculate 
the average value of the modulus of the velocity. To 
leading order in InA^, we find 
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It should be noted that the moduli of the velocity are not 
statistically additive in the sense that 
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This is clear at first sights since the r.h.s. of Eq. (g) 
diverges linearly with the size of the domain. It is 
therefore of interest to calculate the average contribu- 
tion (A|V|(i?)) to the modulus of the velocity exerted 
by vortices located in the annulus of radius R and thick- 
ness Ai?. A similar quantity was calculated by Agekyan 
JlSf in the astrophysical context (see also Kandrup |lj] 
for a generalization of his results to an inhomogenous 
medium) and it is possible to extend his method to the 
case of point vortices. As we shall see, it is possible to 
derive an explicit expression for this quantity in terms of 
Bcssel functions. 

Following Agekyan's method, we proceed in three 
steps: We first determine the average velocity (|V*(i?)|) 
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due to all vortices, except those that lie in the annulus 
AR, assuming that there are N vortices inside R. Then, 
we average this quantity with respect to N. Finally, we 
substract this averaged quantity from the average veloc- 
ity (|V|) due to all vortices. Our presentation closely 
follows the one adopted by Kandrup in his review on 
the stochastic gravitational fluctuations. 

Suppose there are N vortices inside R and A^i between 
R+AR and The probability that these vortices exert 
a velocity V* in O can be obtained by applying Markov's 
method outlined in Ref. |^ . We start from the formula 

Wn,nA^*)= [l[r{r,)d^r,s(v-Y,^X (7) 

i ^ fc ^ 

where T{ri)(Pri governs the probability of occurence of 
the ith point vortex at position r^, and we express the S 
function in terms of its Fourier transform 
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Then, 



with 



W^^.^. = i / ^^.^^ (Py'^^d'p, (9) 
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Let i — 1,...,N label the vortices in the disk of radius 
R (region (I)) and i = N + l,...,Ni label the vortices 
between R + AR and Ri (region (II)). If the vortices are 
uniformly distributed in average, then r(i) = I/ttR^ and 
T(„) = l/7r(i?f - {R^ + Ai?)2). We can write therefore. 
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In the limit Ri, Ni oo with n = Ni/ttR\ fixed (and 
treating InA^i as a constant), we get 



(13) 



with 
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Let P{N) denote the probability of finding TV vortices 
in the disk of radius R. Since the vortices are randomly 
distributed with a uniform distribution in average, P{N) 
is given by the Poisson distribution 
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where n is the average density of vortices. Averaging the 
quantity (|lj) with respect to the Poisson distribution 
(p^), we obtain 



A*(p) = ^P(iV)A^(p) 
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The function C{p) has been calculated in Ref. [^. To 
simplify the integral appearing in Eq. (^6|), we find it 
convenient to take $ as a variable of integration instead 
of r (see Ref. Q for more details). Transforming to po- 
lar coordinates, carrying out the angular integration and 
taking the limit AR — > 0, we obtain 



A*{p) = e-"^('')|l + 2TTnRAR 



, (18) 



where Jq is the Bessel function of zeroth order. Eq. (p^), 
without the subscript N, and Eq. ( [l8| ) determine the 
distribution of the velocity fluctuations occuring in O 
due to all the vortices except those located between R 
and R + AR. If we denote by 



(A|V(i?)|) = (|V|)-(|V* 
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the average velocity due to the vortices located in the 
annulus of radius R and thickness AR, we have 



(A|V(i?)|) 
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Transforming to polar coordinates and carrying out the 
angular integrations, we can rewrite the foregoing expres- 
sion in the form 
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With the change of variables x = 7p/27ri?, we can rewrite 
Eq. (pV) in the form 



(A|V(i?)|) =njAR 
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It turns out that this integral can be expressed in terms 
Under this form, it is not possible to interchange the of modified Bessel functions. Indeed, 
order of integrations. However, using the identity 

d 



xJo{x) = —{xJi{x)), 
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and integrating by parts, we obtain 
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(A|V(i?)|) = 2T:nRAR 
with 

Hp) - pg'ip). (25) 

Integrating by parts one more time with the identity 

4{x) = -J,{x), (26) 

we find 
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For i? — > 0, this expression reduces to 
(A|V(i?)|) ^n-fAR. 



(34) 
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This formula would be expected to be general if the mod- 
uli of the velocity were additive. Indeed, in that case 
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It is now possible to interchange the order of integrations. 
Using the identity 



However, the moduli of the velocity are not additive and 
Eq. (^) does not hold for large values of R. Indeed, for 
R +00, 



Jo{x)dx = 1, 
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we get 
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Straightforward integrations by part lead to the equiva- 
lent expression 



It is seen that "cooperative" effects modify the expression 
of the average velocity that we would expect on the basis 
(29) of naive arguments. In view of the asymptotic behaviors 
( ^ ) and (psi), a simple approximate expression for the 
average velocity is given by 
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Therefore, the average value of the velocity created by 
the point vortices located in the annulus between R and 
R + AR is given by 



l + i?/A' 

where A is a characteristic length defined by 
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This formula shows that the velocity is "shielded" on 
a typical distance A, the inter-vortex distance. There- 
fore, everything happens as if the velocity were additive 
but that each vortex produced an "effective" velocity 



The exact expression for C{p), calculated in Ref. g, is 
relatively complicated. However, for our purposes, it is 
sufficient to consider its approximate expression 



Veff = 
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For R ^ A, the "effective" velocity decays like R ^ in- 
stead of the usual behavior R~^ recovered for i? <C A 
(see Fig. 1). 
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FIG. 1. Plot of the effective velocity normalized by 
the physical velocity V = 'y/2nR versus the normal- 
ized distance R/A. The "effective" velocity is defined by 
Veff = {A\V{R)\) /2nnRAR. The full curve corresponds to 
the exact formula deduced from Eq. ( |34| ) and the dashed line 
to the approximate formula (^l|). 



These results confirms our previous observation 




that the velocity fluctuations occuring in O are domi- 
nated by the contribution of the nearest neighbor. In 
fact, this dominance is marginal since the average veloc- 
ity diverges logarithmically with the size of the domain 

(|V|> - YTR/a'^^ ^ ^n^/\\nR,f/\ (42) 

Therefore, equal logarithmic intervals contribute equaly 
to the velocity. In other words, the contribution from 
vortices between R and -j, -j and and so forth are 
of equal importance. Therefore, the contribution of the 
nearest neighbor (at a typical distance ~ n~^/^ from the 
point under consideration) Vn.n ~ ^v}^"^ is of the same 
order, up to a logarithmic factor, as the contribution of 
the rest of the system. In a sense, we can say that the ve- 
locity ( |4^ ) is dominated by the contribution of the near- 
est neighbor and that collective effects are responsible 
for logarithmic corrections. As described in Ref. ||^, this 
particular circumstance arises because we are just at the 
separation between Gaussian and Levy laws. 

The results obtained in this paper and in Refs. 
can have applications in the context of decaying two- 
dimensional turbulence when the flow becomes domi- 
nated by a large number of coherent vortices. They have 
been used explicitly by Sire & Chavanis to build 
up an effective two- and three-body dynamics of vortices 
subjected to the "noise" induced by the rest of the sys- 



tem. They are also of great importance to build up a ra- 
tional kinetic theory of point vortices in two-dimensional 
hydrodynamics On a formal point of view, they com- 
plete the analogy between 2D vortices and stellar systems 
investigated in Refs. ||l^,|l^. They can also be relevant 
for other fields of physics such as non neutral plasmas 
under a strong magnetic field ||] and parallel straight 
dislocation systems [p^ . 



I acknowledge interesting discussion with C. Sire on 
this topics. 
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